Local theory for ions in binary liquid mixtures 
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Italy 



The influence of ions on the bulk phase behavior of binary liquid mixtures acting as their solvents 
and on the corresponding interfacial structures close to a planar wall is investigated by means 
of density functional theory based on local descriptions of the effective interactions between ions 
and their solvents. The bilinear coupling approximation (BCA), which has been used in numerous 
previous related investigations, is compared with a novel local density approximation (LDA) for 
the ion-solvent interactions. It turns out that within BCA the bulk phase diagrams, the two-point 
correlation functions, and critical adsorption exhibit qualitative features which are not compatible 
with the available experimental data. These discrepancies do not occur within the proposed LDA. 
Further experimental investigations are suggested which assess the reliability of the proposed LDA. 
This approach allows one to obtain a consistent and rather general understanding of the effects of ions 
on solvent properties. From our analysis we infer in particular that there can be an experimentally 
detectable influence of ions on binary liquid mixtures due to steric effects but not due to charge 
effects. 



I. INTRODUCTION 



Ion-solvent mixtures play a central role for various im- 
portant soft matter systems such as colloidal suspensions, 
polymer solutions, biochemical reactions, and electro- 
chemical cells. For most of these systems the presence 
of an appropriate amount of ions is crucial for their func- 
tioning. Therefore the general question arises of how 
and to what extent ions alter the properties of a solvent. 
More than a century ago the considerations by Arrhcnius 
concerning the colligativc properties of ionic solutions [1] 
seemed to imply that the interactions of ions with their 
fluid environment have only a mild effect in the sense that 
adding salt merely modifies the composition of the ion- 
solvent mixture. A few decades later, after more accurate 
measurements had been carried out, this simple picture 
was questioned. Bjerrum, Debye, and Hiickel pointed 
out that electrostatically induced ion-ion correlations are 
expected to influence the ion distribution and the equa- 
tion of state of electrolyte solutions [2, 3]. Since then, 
the Dcbye-Hiickel theory has been widely used, e.g., in 
plasma physics and as an ingredient of the DLVO theory 
(named after Derjaguin, Landau, Verwey, and Overbeek) 
of colloidal suspensions [4] . Within Dcbyc- Hiickel theory 
the solvent is considered to be a uniform dielectric contin- 
uum which influences the Coulomb interaction between 
the ions via a certain permittivity but which is other- 
wise inert. Recently the issue of the mutual influence of 
ions and the solvent has been taken up again, studying 
the solubility of ions and the double-layer structure in a 
near-critical solvent [5-9], possible salt-induced changes 
of the solvent structure [10-17], and effects of the inho- 
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mogencities of the permittivity close to interfaces [18-21]. 

These investigations require a model for the solvent at 
least on the mesoscopic scale as well as a description of 
the ion-solvent interaction. To this end one can split the 
pair potential between the species into the long-ranged 
electrostatic monopole-monopole contribution and into 
the remaining contributions of shorter range, which we 
shall refer to as chemical contributions. In the vast ma- 
jority of the theoretical studies of dilute electrolyte so- 
lutions, ions are described as point-like particles whose 
chemical contributions to the interactions with the sol- 
vent are modeled locally within the so-called bilinear cou- 
pling approximation (BCA). This amounts to a local den- 
sity approximation for the chemical contribution to the 
excess free energy which is bilinear in the particle number 
densities [5 8. 10 12, 18, 19, 21]. Within the approaches 
of Refs. [9] and [20] the ion size is accounted for by means 
of hard-core exclusion and solvation is modeled by non- 
local interactions within the so-called random-phase ap- 
proximation (RPA) for density functional theory (DFT). 
In fact, the BCA can be considered as the local version of 
the RPA, which is expected to be reliable only for inter- 
action energies small compared with the thermal energy 
[22]. However, the ion-solvent interaction is typically of 
the order of some tens of the thermal energy [23, 24]. 
Therefore the application of the BCA or of the RPA to 
ion-solvent mixtures is questionable [25]. 

In the following we show that for realistic values of the 
parameters the BCA and the RPA indeed lead to unphys- 
ical results. Since non-local models are notoriously com- 
plicated it seems worthwhile to investigate the possibility 
of local descriptions of ion-solvent interactions. In order 
to demonstrate that valuable improvements within the 
class of local models are possible, we propose an alterna- 
tive local density approximation (LDA) the predictions of 
which are in qualitative agreement with experimental re- 
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suits and which does not lead to the artifacts introduced 
by the BCA. Without presenting the details which will 
be expounded below, this LDA has already been applied 
successfully in a recent study of the effective interaction 
between two planar substrates in contact with a near- 
critical binary liquid mixture and in the presence of salt 
[26]. As far as the phase behavior, the bulk structure, 
and the asymptotic intcrfacial structure of the solvent 
are concerned, our results are equivalent to replacing the 
bare interaction energies used within the BCA by effec- 
tive ones which saturate for large values of the bare inter- 
action energies. Within this LDA, and in contrast to the 
BCA, more realistic estimates of the magnitude of salt- 
induced effects can be obtained, which are expected to 
be important for interpreting and designing experiments 
and for considering applications. 

In the following, after introducing the model and the 
LDA in Sec. II and in Appendix A, bulk systems are dis- 
cussed in Sec. Ill, where we focus on the phase diagram 
and on the structure of the correlation functions. Similar- 
ities and differences between the proposed LDA and the 
BCA are highlighted and compared with the available ex- 
perimental data. In Sec. IV intcrfacial structures as well 
as critical adsorption in semi-infinite planar systems are 
discussed. For critical adsorption important qualitative 
differences between the LDA and the BCA are revealed, 
with the predictions of the former being in agreement 
with the corresponding qualitative experimental findings, 
whereas those of the latter are not. We propose further 
experimental investigations which are expected to dis- 
criminate more sharply between the LDA and the BCA 
than the presently available data do. In Sec. V an appli- 
cation of the proposed LDA to colloidal interactions in 
near-critical electrolyte solutions is discussed and com- 
pared with alternative approaches described in the liter- 
ature. Finally, in Sec. VI we draw our conclusions and 
provide a summary. 

II. MODEL 

A. Definition 

We consider a three-dimensional (d = 3) container 
VCR 3 filled with an incompressible binary liquid mix- 
ture acting as a solvent for cations (+) and anions (— ). 
All solvent particles are assumed to be of equal size with 
non-vanishing volume a 3 whereas the ions are consid- 
ered to be point-like; hence ions do not contribute to 
the total packing fraction (see also Appendix A). The 
set of dimensionless positions r = (x,y,z) := r/a for 
r G V is defined as V. At r G V the number densi- 
ties of the solvent components A and B are given by 
qa(^) = 0(r)a~ 3 and Qb(t) = (1 — </>(r))a~ 3 , respec- 
tively, with < </> < 1, whereas the number densities of 
the cations and anions are given by g + (r) = g + (r)a~ 3 
and Q-{v) = g^(r)a~ 3 , respectively The walls dV of 
the container carry a surface charge density a(r)ea~ 2 at 



r G dV, where e is the (positive) elementary charge. The 
influence of the walls onto the solvent due to short-ranged 
chemical effects is captured by surface fields localized at 
the walls. At r G dV, the dimensionless volume fraction 
4>(r) of A particles couples linearly to surface fields h(r), 
where h > (< 0) leads to a preferential adsorption of 
solvent component A (B). The equilibrium profiles <fi, 
g + , and g_ minimize the approximate grand potential 
density functional ksTQ[4>, g±], 

{![</>, Q*] = J^d 3 r jawMr)) + ^W(r)) 2 

+EH™ (ft(r)) + ft(r)KiWr)) ] (1) 

i=± 

+ ^T D ( r '[^) 2 i- / d ^Mr)0(r), 

withuw^) = 0(]n0-^) + (l-0)]n(l-^) + x(T)0(l- 

</>) and u\^{g±) = f?±(ln g± — 1 — /j,±) as the bulk grand 
potential densities of the solvent and of the ±-ions (in 
the low number density limit), respectively. Here ksT is 
the thermal energy, fi^ksT and ^,±ksT are the chemical 
potential difference (fiA — ^s)fcsT of the solvent par- 
ticles and the chemical potentials of the ±-ions, respec- 
tively, and Ibcl = e 2 / (47r£ofcsT) is the Bjerrum length for 
the vacuum permittivity sq. The temperature-dependent 
Flory-Huggins parameter x(T) > describes the effective 
interaction between solvent particles, where the temper- 
ature dependence is usually described by the empirical 

Xh 

form x(T) = xs H — with the system specific entropic 

contribution xs an d the cnthalpic contribution xh [27]. 
For x{T) > x{Tc) phase separation occurs in the pure, 
salt-free solvent within a certain range of <p whereas for 
x(T) < x(^c) the solvent components A and B are mis- 
cible in any proportion. A positive (negative) enthalpic 
contribution xh corresponds to an upper (lower) crit- 
ical demixing point. The gradient term cx (V</>(z)) 2 
with V = aV penalizes the spatial variation of the sol- 
vent composition [28]. The ion-solvent interaction is de- 
scribed within a local density approximation (LDA) by 
the effective ion potential fcsTV±(^>) generated by the 
solvent (see below). The relative permittivity e((j)(r)) 
is assumed to depend locally on the composition of the 
solvent (j)(r) but not on the ion densities g±(r), which 
is justified for small ionic strengths, i.e., g±(r) <C 1. 
Here the mixing formula s(4>) = ea4> + £ b(1 — 0) intro- 
duced by Bottcher [29] is used [5, 19, 21]. Using Si-units, 
the electric displacement D = Dea -2 in Eq. (1) fulfills 
Gauss' law V ■ D(r, [g±]) = Q+(r) — g_(r),r G V, with 
fixed surface charges n(r) ■ D(r, [g±]) — er(r),r G dV, 
where n is the unit vector perpendicular to dV point- 
ing towards the exterior of V (see Ref. [4]). Note that 
D(r, [g±]) is generated by the ±-ions and the given sur- 
face charges a; it does not depend explicitly on <fi. Within 
the present model, besides being confined, ions interact 
with the walls only electrostatically. 

Note that by using the square-gradient form of Eq. (1) 
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FIG. 1: Bulk phase diagrams of a binary liquid mixture with added salt of constant chemical potential (per ksT) hi = /i+ + /i_ 
within the bilinear coupling approximation (BCA) V± C (4>) = f±4> [( a ) an d (b)] and within the local density approximation 
(LDA, see Eq. (3)) V±(<f>) = — ln(l — <j>(l — exp(— /±))) [(c) and (d)] in terms of the Flory-Huggins parameter \ and the 
composition (f> [(a) and (c)] or the chemical potential (per fcsT) /x^ = ha — fJ-B conjugate to the composition cf> of the binary 
solvent [(b) and (d)]. The thick solid lines correspond to the binodals, which delimit the two-phase coexistence regions in 
the 4>~x diagrams [(a) and (c)] from below. The dashed lines are the spinodals and the thin horizontal line in panel (a) 
is the tie-line corresponding to the triple point (a) found within BCA. Representative values for the solubility contrasts 
per ksT, (/+,/-) = (3,26), have been chosen. The chemical potential \iiksT of the salt corresponds to an ionic strength 
I c = Q± — 10 mM at the critical point with composition 4> = 4>c ~ § • The weak influence of the salt on the phase diagram 
within LDA leads to curves in panels (c) and (d) which are, on the present scale, almost (but not quite) symmetric with respect 
to <f> = | and H4> = 0, respectively. Whereas the LDA [(c) and (d)], in agreement with the experimental evidence, exhibits a 
single critical point (•, </> c> i « ~, Xc,i ~ 2), which slightly shifts upon changing the ionic strength (see Fig. 3), the standard BCA 
[(a) and (b)], in contrast to the available experimental observations, leads to a second critical point (•, cj>c,2 ~ 0.1, Xc, 2 ~ 2.1) 
as well as to a triple point (A). 



we implicitly assume that the interactions are short- 
ranged [30], i.e., van der Waals forces are not taken into 
account. Moreover, layering due to packing effects close 
to walls is also not accounted for by square-gradient the- 
ories. Nonetheless such a description provides reliable re- 
sults at mesoscopic scales [30] . Finally, the ionic strength 
is assumed to be sufficiently low so that one can neglect 
short-ranged ion-ion interactions. Therefore the ions in- 
teract with each other only via the electrostatic field. 



Accordingly, the expression for w ion does not contain 
additional Flory-Huggins parameters and there are no 
square-gradient terms for g±. However, these features of 
the simple functional in Eq. (1) are not expected to lessen 
the main conclusions of our study, which is devoted to 
investigate the kind of influences of ions on solvent prop- 
erties, rather than to construct models with quantitative 
predictive power. 
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B. Ion-solvent interaction 

In Eq. (1) the ion-solvent interaction is described, 
within a local density approximation (LDA), by a 
solvent-induced ion potential, V±((/>)fcsT. The bilinear 
coupling approximation (BCA) used in previous inves- 
tigations (see, e.g., Refs. [5-8, 10-12, 19, 21]) corre- 
sponds to the choice V^ CA \(j)) :— f±<p, where f±ksT = 
(f±A — f±B)ksT is the difference between the bulk solva- 
tion free energies of a ±-ion in solvents consisting purely 
of component A, f±AksT, and purely of component 
B, f±BksT. The solubility contrasts f±ksT are also 
known as Gibbs free energies of transfer. In this con- 
text the only relevant parameters arc the two differences 
f± = f±A~f±B because the other two independent quan- 
tities f±A~\~f±B can be absorbed as shifts in the definition 
of the chemical potentials fi±ksT of the ions. For bulk 
systems the BCA is identical to the random phase ap- 
proximation (RPA) [22] , which is expected to be reliable 
only if the coupling strengths are much smaller than the 
thermal energy, i.e., \f±\ <C 1. However, for electrolyte 
solutions, this condition is in general not fulfilled. In- 
stead, the Gibbs free energies of transfer between two 
liquids are usually of the order of some lOfcgT [23, 24]. 

Figures 1(a) and 1(b) display the bulk phase dia- 
gram for a constant chemical potential (per ksT) \ii := 
/i+ + /i_ of added salt obtained within BCA for the rep- 
resentative values /+ = 3, /_ = 26. This choice is similar 
to the Gibbs free energies of transfer for potassium chlo- 
ride (KC1) from water to acetone: /+ = 2, /_ = 23 [23]. 
The condition of local charge neutrality Q+ = Q- =: I in 
the bulk implies that the ionic strength I depends on the 
chemical potentials /i±fcsT of the ions only via the sum 
= /!/. For given uniform composition </> and ionic 
chemical potential \ii the Eulcr-Lagrangc equation of 
in Eq. (1) with respect to uniform ion densities g± = I 
can be used to express the bulk ionic strength as 
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Ibul] 



b, fij) = exp - V+W) - V- m) . (2) 



Within the present model, /bulk is independent of the 
Flory-Huggins parameter Xi i- e -i h depends on the tem- 
perature T only via the normalizations of \ii and V±, 
which arc defined in units of ksT. In Fig. 1 the chem- 
ical potential [ii of the salt is fixed such that the sol- 
vent composition (j> — \ leads to an ionic strength 
/ = Ia~ 3 = 10 mM ps 0.006 nm -3 , where here and in 
the following we choose the length scale a — 2 A. Due to 
the absence of gradients, electric fields, and surfaces, the 
bulk phase diagram is determined by the first, third, and 
fourth term on the right-hand side of Eq. (1). The the- 
oretically predicted occurrence of two critical points (•) 
as well as of a triple point (A) is not supported by exper- 
imental evidence, which signals the breakdown of BCA 
for such large parameters f± . Whereas for most systems 
it is experimentally difficult to preclude the occurrence 
of such a second critical point or triple point, the experi- 
mental resolution is yet sufficiently high to exclude these 
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FIG. 2: Comparison of the solvent-induced ion potential 
V±(4>) [(a)] and its derivative j± = V±((j>) [(b)] within LDA 
and BCA for ion solubility contrast f± (see the main text). 
For small values of f± (see the case f± = 1) the differences 
between LDA and BCA are small. For large values of f± (see 
the case f± = 10) V± (<j>) and 7± = V± (4>) become large at sol- 
vent compositions (f> ps 0.5 within BCA whereas they remain 
small within LDA. Within LDA V±((f> = 0) = 1 - exp(-/±) 
and V±(4> = 1) = exp(/±) - 1, while within BCA V±(<£) = f±. 



features to occur visibly to the extent as predicted by the 
BCA (Figs. 1(a) and (b)). 

A more appropriate approximation for the solvent- 
induced ion potential V±(0)fcsT, which is derived in Ap- 
pendix A, is given by (see also Fig. 2(a)) 



V± (</>) = - ln(l - 0(1 - exp(-/±))). 



(3) 



F° r l/±l *C 1 this expression reduces to the correct 
asymptotic expression V± {<p) ~ t/_( bca ) ^ _ j n nrm t 
f± — > oo, i.e., if ions are insoluble in component A, 
V±{4>) ~ — ln(l — 0), which corresponds to the free en- 
ergy of the ions dissolving entirely in component B only, 
which has the volume fraction 1 — 0. Similarly, in the 
limit f± — > — oo, i.e., if the ions are insoluble in compo- 
nent B, V±(4>) ~ f± — ln0, which is the free energy of 
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the ions dissolving entirely in component A only, which 
has the volume fraction and for which the solvation 
free energy is f±. For the same set of parameters as in 
Figs. 1(a) and (b), Figs. 1(c) and (d) display the phase 
diagram within LDA. In agreement with experimental 
observations, within LDA only a single critical point (•) 
occurs (see, e.g., the closed loop-binodals in Rcf. [17] with 
only one lower critical demixing point in the presence of 
an antagonistic salt, i.e., with /+ and /_ having opposite 
signs). Hence one can conclude that the standard BCA, 
i.e., V( BCA \cj)) = f±(/), introduces artifacts for too large 
ion-solvent couplings, 3> 1, which are absent within 
the LDA proposed in Eq. (3). Note that within LDA the 
salt has such a weak influence on the phase diagram that 
the curves in Figs. 1(c) and (d) are, on that scale, almost 
(but not quite) symmetric with respect to = , an( i 
/j0 = 0, respectively, whereas within BCA the influence 
of salt on the phase diagram is so strong that the curves 
in Figs. 1(a) and (b) are highly asymmetric. 

If deviates slightly from a certain composition 
0o G [0,1] one has V±(0) ~ V±(<fo) + 7±(0 - 0o) 
with the effective coupling strengths j± := V±(<fro) — 
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1 - exp(-/±) 



1 



1 



instead of 



1 - 0o (1 - exp(-/±)) L o 'l-0 O - 
f± as in BCA. For, e.g., 0o = 1/2 one finds j± — 
2tanh(/±/2) € [-2,2], i.e., the use of BCA, which cor- 
responds to 7± w f±, is justified only for Gibbs free 
energies of transfer per fc^T, f±, not larger than 2 (see 
Fig. 2(b)). However, in previous investigations BCA has 
been used even for large values of |/±| [5-8, 19, 21]. 
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III. BULK SYSTEMS 

Bulk properties such as the phase diagram and par- 
tial structure factors are determined routinely in order to 
characterize experimentally the behavior of fluids. The 
available experimental data offer the possibility to as- 
sess the quality of the proposed LDA (see the previous 
Sec. II) with respect to predictions of bulk properties and 
to compare it with the frequently used BCA. 



FIG. 3: Variation of the critical volume fraction C [(a)] and 
the critical Flory-Huggins parameter Xc [(b)] as function of 
the ionic strength I c = I c a~ 3 at the critical point for two rep- 
resentative sets of solubility contrasts: (/+,/_) = (3,26) and 
(0, 20). These results show that both | — C and \c — 2 depend 
linearly on I c and that there are no quantitatively significant 
shifts of the critical point upon varying the ionic strength 
within experimentally reasonable ranges. On this scale, the 
phase diagrams for (/+, /_) = (3, 26) (see Figs. 1(c) and (d)) 
and for (/+,/-) = (0,20) are almost indistinguishable. Note 
that C (J C -»• 0) = i and X c{L -» 0) = 2. 



A. Critical point 



There is experimental evidence [31] that in the phase 
diagram of a binary liquid mixture the critical point 
shifts upon adding salt. The direction as well as the 
magnitude of the shift depend on the materials prop- 
erties of the binary liquid mixture and of the ions. 
Due to the relations qa = 0, Qb = 1 — 0, Q+ = 
Q- = I the bulk system, which comprises four par- 
ticle species, is de facto a binary mixture, character- 
ized by {10, \ij = /U+ + /j,-, and x (i- e -> T). Hence 
in this three-dimensional space of thermodynamic vari- 
ables there is a sheet of first-order demixing phase tran- 
sitions (fJ^ sm (l^i, x), M/j x) bounded by a line of crit- 



ical points (/i^ crlt ^(/zj),/i7, X^ cllt H^i)) which translates 

into (0 c ,/ c ,Xc) = (0 (crit) (/i/),/ (crit HM/),X (cr " ) (M/))- For 
a given chemical potential fiiksT the critical point 
(0o Ic, Xc) is determined as the minimum of the Flory- 
Huggins parameter Xs(0> -^buik(0, A*/)) a t the spinodal as 
a function of for constant /i/. The spinodal is defined 
by the set of points (0, /, Xs(0, 1)) m the bulk phase dia- 
gram for which points (0, /, x) with x > Xs(0i I) exhibit 
no longer at least a local minimum of the density func- 
tional Eq. (1) (see the dashed lines in Figs. 1(a) and (c)). 
Accordingly, at the spinodal the Hessian matrix of the 
bulk grand potential density 0(0, 1)/V corresponding to 
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Eq. (1) has a zero eigenvalue. This condition leads to 

1 



1 



(4) 



n 7 w + vi (0) 



By inverting the relation I c = /( crlt ) (^tj) one obtains /ij = 
/ij C1It '(/ c ). Figure 3 displays the variation of (a) the crit- 
ical volume fraction 4> c = <^ crlt ' (^j (ic)) and (b) the 
critical Flory-Huggins parameter \ c = X^^i^i (Ic)) 
as functions of the ionic strength I c = I c a~ 3 at the criti- 
cal point within the present LDA model. Without added 
salt (J = 0) one obtains the critical point (<f) c , Xc) = (|, 2) 
of the pure solvent. For the given choice of the parame- 
ters (/+, /_) and for small I c the critical composition <p c 
decreases and the critical Flory-Huggins parameter \c in- 
creases linearly upon increasing the ionic strength I c . For 
small I c the asymptotically linear dependence of the crit- 
ical point (0 C , Xc) on the ionic strength I c is in agreement 
with experimental evidence [31, 32]. However, the mag- 
nitudes of these shifts are tiny, even for large differences 
in the solubility contrasts, e.g., (/ + , /_) = (0, 20), to the 
effect that the bulk phase diagrams are almost indistin- 
guishable within experimentally relevant ranges of ionic 
strengths I c < 10 mM. Although for the parameters used 
in Fig. 3 the shifts of the critical point (4> C iXc) as pre- 
dicted within the BCA are up to three orders of magni- 
tude larger than within the LDA, the effect is still small. 
Within the range of ionic strengths considered in Fig. 3 
the experimentally observed critical point shifts are also 



small [31]. However, significant shifts of the critical tem- 
perature T c have been detected for large ionic strengths 
I c ^s> 100 mM [31]. But for such large ionic strengths the 
model in Eq. (1) is not expected to be applicable, because 
it neglects short-ranged ion-ion interactions. 



B. Correlation functions 

The bulk structure of fluids, which is experimen- 
tally accessible by X-ray and neutron scattering, pro- 
vides information complementary to those which follow 
from the bulk phase behavior. Hence it provides addi- 
tional opportunities to assess the quality of the LDA. 
Here we consider a spatially uniform equilibrium state 
{(f), I, x) m the one-phase region of the phase diagram 
(see Fig. 1(c)), which minimizes the density functional 
f2 in Eq. (1) in the absence of surfaces, i.e., without 
the last term therein. The corresponding two-point cor- 
relation functions Gij{r) = QiQjhij(r) = QiQj(gij(r) — 
G {</>, +,— },Q4> ■= (f>, Q± = I, are obtained from 
Gij(r) = Gij(r,0),r = |r|, with the inverse G~Mr, r') = 



where 



E 

3 



dV 



G^(r,r')G jk (r',r") 



SQi(T)Sg j {r , y — Jv 
5ik5(r — r ). The three-dimensional Fourier transforms 

oo 

4-7T 

T 



Gij(k) 



dr rGij (r) sin(fcr) with dimensionless k, 



which arc proportional to the partial structure factors 
[22], are given by 



J 



G^(fc) 
G 0± (fc) 
G±±(k) 
G± T (k) 



1 



L(k) 
I 



(k 2 + K 2 



L(k) 
I 



1 



1 - 
1 



2 X +i{v"{4>)+vnm 



L(k) \\(/) 1 - 



+ f^ 2 ) ^ + Vi(cf>)V^) 2 Ik 2 



K 

~2~ 



V' ((f>) 2 Ik 2 



K 

T 



(5) 



with Note that V±(4>) = leads to G^±(k) = 0, i.e., as ex- 

8w£bI pected, the fluctuations of the solvent composition and 

(6) of the ion densities are uncorrclatcd in the absence of 
ion-solvent interactions. 



K 2 



e(<f>) 



as the square of the inverse Debye length and the denom- 
inator (see Eq. (5)) Due to the constraint qa + Qb = 1, the correlation 
/ X \ functions Gaa{^), G J 4s(fc), and Gss(fc) of the number 
L(k) := (k 2 + k 2 ) i^ — k 2 + 2(xs(4>',I) ~ x)J density fluctuations of the A and B particles are re- 
j lated to the correlation function G^{k) by GAA(k) = 
--(V+(</>) - V_{(f>)) 2 k 2 . (7) ~G AB {k) = G BB {k) = G^k). In the following we re- 
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fer to G^(k) as the solvent structure factor. It can be 
written in the form 



G^(fc) 



G W(0) 



1 + (Xk) 1 - 



.9 



1 + (*/«) 



(8) 



with 



and 



6(x s (<£,/)-x) 



2 3(A 7 ) 2 eW 
.9 := 



(9) 



(10) 



where A7 := 7+ - 7- = V\_((f>) - V!_{4>). The isother- 
mal compressibility, which is proportional to G^(0) = 
(2(x s (0, /)-x)) _1 , diverges oc \x~Xc\~ 1 upon approach- 
ing the critical point (<f> c , I Cl x c )- As expected within the 
present mean-field theory, one finds the classical criti- 
cal exponent 7 = 1 instead of 7 as 1.24 for the Ising 
universality class [33] . For a state point ((f), /, \) m the 
bulk phase diagram (see Figs. 1(a) and (c)) the length 
A is an (inverse) measure of the deviation of x from its 
value Xsifiil) at the spinodal. Equation (8) has already 
been derived in Ref. [5] within BCA, which corresponds 
to the linear approximation A7 sa /+ — /_. For \g\ < 1 in 
Eq. (8) the solvent structure factor G^(k) is a monoton- 
ically decreasing function of the wave number k, whereas 
for \g\ > 1 at k 
G^(0) 



= ^Vu/j-l a maximum G^(fc max ) = 
2 - occurs. Hence, if \g\ > 1, G^(k max ) 



\-(nX) 2 (\g\-l) 2 
diverges as function of x at A = A unst abio = - 1)) \ 

i.e., the spatially uniform bulk state becomes unstable 
upon approaching the critical point. Note that in the 
limits \g\ — > (no ion-solvent coupling) or n — > (no 
salt) Eq. (8) leads to the Ornstein-Zcrnikc-like solvent 
structure factor G^(k) = G<^(0)/(1 + (Xk) 2 ). In this 
case A can be identified with the bulk correlation length. 

Experimental reports of uniform bulk states close to 
the critical point of water+2,6-dimcthylpyridinc mix- 
tures with KBr, KC1, and Mg(N0 3 ) 2 (see Ref. [34]) as 
well as distributions of neutron scattering intensities of 
water+3-methylpyridine with LiCl, NaCl, KC1, NaBr, 
and MgSC>4, which vary monotonically as function of k 
(see Refs. [13, 14]), indicate that in these systems one 
has \g\ < 1. Within the present LDA this latter rela- 
tion is expected to be fulfilled: Close to the critical point 
(0oXc) ~ (|,2) (see Fig. 3) of, e.g., the widely studied 
binary liquid mixture of 3-methylpyridine (component A, 
ea = 10) and water (component B, sb = 80) with a lower 
critical demixing point at T c rj 316 K, i.e., ps 529 A, 
one obtains \g\ < 0.3 independent of the type of salt, be- 
cause |A 7 | < 4 (see the last paragraph in Subsec. IIB). 
However, within BCA, i.e., for A7 xs f + — f_ with typ- 
ically |/ + - /_| > 1 [23, 24], one has to expect \g\ > 1, 



which, according to the above reasoning, is in sharp con- 
trast to the available experimental results. We men- 
tion that experimental reports [15-17] of "periodic struc- 
tures" in heavy water+3-methylpyridine mixtures with 
sodium tctraphenylborate (NaBPh4) cannot, however, be 
expected to find a consistent interpretation in terms of 
a local ion solvation model as given in Eq. (1), neither 
within BCA nor within the present LDA, because the 
anions ([BPli4] _ ) are much larger than the solvent par- 
ticles, such that in these systems the ion size is expected 
to be relevant. 

The charge-charge structure factor Szz(k) = 
(G±±(k) — G±if(k))/I [22], which measures correlations 
of fluctuations Z of the local charge density around 
g + — p_ = 0, is obtained by inserting the expressions 
for G±±(k) and G± T (k) from Eq. (5): 



Szz(k) = k 



2 3 



k 2 + 2( X s(^I)-x) 



L(k) 



(11) 



The asymptotic behavior Szz(k — > 0) ~ (k/n) 2 is the 
signature for perfect screening [22]. Further, the case 
A7 = corresponds to the Dcbyc-Hiickel limit Szz (k) = 

k 2 /(k 2 + K 2 ). 

The asymptotic behavior of the correlation function 



Gij(r) 



1 



2ir 2 r 



dk kGij(k) sin(fcr) can be inferred from 



a pole analysis of G,j (k), which amounts to determine the 
roots of the denominator L(k) defined in Eq. (7) [35, 36]. 
Since L(k) is a polynomial in k of degree four it has 
four and only four complex roots fc„ = k' v + ik",k' v = 
Re(£v),/c" = ~kn(k u ),v G {1,...,4}. Due to the actual 
structure of L(k) there are constraints on the locations 
of the four roots fc„ in the complex plane. If L(k = k v ) 
vanishes this holds also for k = k*, because L(k) has 
real coefficients. Moreover, if L{k = k u ) vanishes this 
also holds for k = —k Vl because L(k) is a polynomial 
in k 2 . Accordingly this is also true for k = — k*. This 
implies the root structure shown in Fig. 4. Three dis- 
tinct situations can occur. For purely imaginary roots 
given by {fci = ik",k2 = ik'^k^ = —k\,k^ — —k^ with 
< k'l < k'2 (see Fig. 4(a)) the asymptotic decay of the 
two-point correlation functions Gjj (r — > 00) is monotonic 
(x exp(— k'{r)/r. For complex roots {k\ = k[ + ik", &2 = 
-kf,k 3 = -ki,ki = k{\ with k[,k'{ > (see Fig. 4(b)) 
the two-point correlation functions Gij(r) vary asymp- 
totically oc sm(k[r + const) exp(—fc"r)/r giving rise to 
a damped oscillatory decay. Finally, purely real roots 
{k\ = k'i, = k' 2 , &3 = —ki, k± = —k 2 } with < k[ < k' 2 
(see Fig. 4(c)) indicate an unstable bulk state, i.e., the 
corresponding point in the phase diagram is located in 
between the spinodals. The exponential decay of the 
two-point correlation functions (whether monotonically 
or oscillatory) is consistent with the short range of the 
interactions implied by taking a gradient expansion in 
Eq. (1). 
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FIG. 4: Poles fci, . . . , &4 of the Fourier transform Gij (fc) of the 
two-point correlation functions Gy(r) in the complex plane 
k = k' + ik" G C, which correspond to the roots of the denom- 
inator L(k) (see Eq. (7)). According to the analytic structure 
of L(k) (see the main text) only the three distinct situations 
shown in panels (a)-(c) can occur. Purely imaginary poles 
[(a)] correspond to a monotonic decay of Gij(r — > oo) whereas 
a pole structure as in panel (b) (|fc„| all equal) corresponds 
to an oscillatory decay of Gij(r — > oo). Purely real poles 
[(c)] indicate an unstable bulk state, which does not occur in 
the one-phase region of the phase diagram in Fig. 1(c). The 
merging of two poles on the imaginary axis [(d)] corresponds 
to a Kirkwood crossover point. 



Thermodynamic states in the bulk phase diagram with 



monotonically decaying Gi 



oo) are separated from 



states with damped oscillatory decay of Gy(r — > oo) by 
so-called Kirkwood crossover lines [37]. Crossing these 
lines is associated with the merging of two purely imagi- 
nary poles (see Figs. 4(a) and (d)) of Gij(k) in the upper 
(and similarly in the lower) half of the complex plane 
and with a subsequent emergence of a pair of two poles 
(see Fig. 4(b)) with equal imaginary parts and with real 
parts of equal absolute value but of opposite sign [38]. 
In the phase diagrams of Fig. 5 the Kirkwood crossover 
lines are denoted by dotted lines and damped oscilla- 
tory decay of Gij(r — > oo) occurs at state points in the 
grey area enclosed by the Kirkwood crossover lines. The 
parameters correspond to the aforementioned binary liq- 
uid mixture water+3-methylpyridine (ea = 10, Sb = 80, 
£sa = 529 A); for simplicity the temperature depen- 
dence of the Bjerrum length £g is ignored. The chem- 
ical potential [ii of the salt is fixed such that there is 
an ionic strength I c = 10 mM at the (shifted) critical 
point. Figures 5(a) and (c) correspond to the parame- 
ters (/+,/_) = (3,26) used in Figs. 1(a) (BCA) and (c) 
(LDA), respectively. Figure 5(b) refers to the case of a 
strongly antagonistic salt, (/+,/-) = (—20,20), whereas 



Fig. 5(d) relates to the intermediate case (/+,/_) = 
(0,20). Within BCA (see Fig. 5(a)), the damped oscilla- 
tory decay of Gy(r — > oo) prevails in a large portion of 
the phase diagram, and wave lengths of the oscillations 
as small as the particle size a can occur at state points in 
the center of the grey area. However, within the present 
LDA, damped oscillatory decay of GV, (r — > oo) is found 
only in a narrow range of O(A'yl) for values of x around 

_ Xs{<P,I) o((A"f) 2 I), which extends into the one- 
1 + k 2 /6 u ' ' 
phase region only in the vicinity of the critical point (see 
Figs. 5(b)-(d)). 

For small wave numbers k the structure factor 
G<^(0) 



G^ik < k) 



Ornstein-Zcrnikc form 



l + (AA0 2 (l- 9 2 ) 
G</></>(0) 



l + (£(OZ) fc ) 



(see Eq. (8)) takes the 
with the length 



I x(i-g 2 ) 

6( Xs (^I)-x) 



(12) 



which we shall refer to as the Ornstein-Zernike length. 
This length £( oz ) is determined routinely in scatter- 
ing experiments by fitting an Ornstein-Zernike expres- 
sion to scattered intensities at small momentum trans- 
fer [34]. For water+2,6-dimethylpyridine mixtures with 
KBr, KC1, and Mg(N0 3 ) 2 (see Ref. [34]) it has been 
found experimentally that the amplitude of £( oz ) = 

£o K-^ 1 — T c )/T c \~ v is to a large extent independent of 
the considered type of salt and ionic strength. Due to 
^q 02 ' oc yjl — g 2 this observation indicates that one has 
g 2 <C 1, which, according to the arguments given above, 
is expected within LDA but is not compatible with pre- 
dictions following from BCA. 

The poles {ki, . . . , k^} of the solvent structure factor 
G ( f >( f > (k) can be expressed in terms of the Ornstein-Zcrnikc 
length £( oz ) and the inverse Debye length k. For a mono- 
tonic decay of G^^r — > oo) one has purely imaginary 
poles k v = ik'l with 



k'{ = -4' = \j2(u- ^~^) 



and 



k'{ = -k'i = ^2(u + v^ 7 ^), 



(13) 



(14) 



whereas a damped oscillatory decay of G<p<p{r — > oo) 
is characterized by the poles at {k\ = k[ + ik",k 2 = 
— fcjf, fc 3 = —k±, k± = k*} with 



k[ = \/v — u, k'l = \/v + u, 



(15) 



where 



n 2 (l-g 2 ) 



( 1 + R<M (16) 
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FIG. 5: Phase diagrams as in Fig. 1 with Kirkwood crossover lines (dotted lines) within the bilinear coupling approximation 
(BCA) [(a)] and the local density approximation (LDA) [(b)-(d)]. The parameters correspond to the binary liquid mixture 
water+3-methylpyridine (ea = 10, eb = 80, Ibcl = 529 A); for simplicity the temperature dependence of the Bjerrum length £b 
is ignored. The chemical potential fxiknT of the salt is fixed such that at the (slightly shifted) critical point with composition 
(f> = (f> c « ~ there is an ionic strength I c — 10 mM. Outside the grey regions bounded by the dotted lines the two-point 
correlation functions exhibit asymptotically a monotonic decay, whereas inside these regions damped oscillatory decays occur. 
Within BCA a large portion of the phase diagram corresponds to oscillatory decay, whereas within LDA this occurs only in 
a narrow band within which the value of the bulk correlation length is close to that of the Debye screening length. Note the 
differences in scales for the axes in (a) and in (b)-(d). 



and 



2£(° z ) ' 



(17) 



Close to the critical point (i.e., for £( oz ) — > oo) one finds 
a monotonic decay of G^(r — > oo) with the decay length 
l/k'{ ~ £< oz ) oc \x~Xc\~ v (see Eq. (13)) with the mean- 
field critical exponent v =■ \ instead of v rj 0.63 for the 
Ising universality class [33]. Therefore the electrostatic 
interactions do not affect the universal critical exponent, 
but they can influence the non-universal critical ampli- 
tudc ^° Z) . 

Figure 6 displays the real and imaginary parts of the 
poles k v of Gij(k) in the ranges k' v = Re(/c J/ ),fc" = 



Im(fc„) > at the critical composition <fi = <fi c « i for the 
parameters corresponding to Fig. 5(b). The four poles 
{ki, . . . , £4} can be expressed in terms of (k[, k'{, k'2) (see 
Eqs. (13)— (15)). In Fig. 6 the two purely imaginary 
poles ki and ki with positive imaginary parts for mono- 
tonically decaying Gy (r — > 00) occur as two branches, 
which merge at the Kirkwood crossover points (•). From 
Eqs. (13) and (14) the Kirkwood crossover points (•) are 
characterized by u = v which leads to (see Eqs. (16) and 



< 17 » < (oz, ^ 



1«1 Hence at the Kirkwood 



crossover points the inverse decay lengths of Gij(r — > 00) 
correspond approximately to k (dashed line) and l/^ oz ' 
(dotted line). At the critical point (x = Xc, i-e-, k" = 0, 
see Eqs. (13) and (17)), Gy(r — » 00) decays as 1/r with 
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FIG. 6: Within LDA real and imaginary parts of the poles 
k v = k' v + ik",k'„ = Re(/c„),fc" = Im(fc„),f € {1,...,4}, 
of the Fourier transform Gij(k) of the two-point correlation 
functions Gij(r) as functions of the deviation Xc — X from 
the critical point at the critical composition cf> = (f> c w | 
for the parameters corresponding to Fig. 5(b). The four 
poles ki, . . . , ki can be expressed in terms of (k[, k" , k'2) (see 
Eqs. (13)-(17) and Fig. 4). If Gij(r — > 00) decays mono- 
tonically, the poles of Gij(k) are purely imaginary (k[ = 0), 
giving rise to two branches k'{ and k'2 of positive imaginary 
parts (see Fig. 4(a)). If Gij(r — > 00) decays oscillatorily, 
there is only one pole of Gij(k) with positive real and imagi- 
nary parts (k'2 = k") (see Fig. 4(b)). The merging of the two 
branches k'{ and k'2 for monotonic asymptotic decay takes 
place at the Kirkwood crossover points (•) (see Fig. 4(d)). 
Upon varying <f> these points form the Kirkwood crossover 
lines (dotted lines in Fig. 5). For comparison the inverse De- 
bye length k (dashed line) as well as the inverse Ornstein- 
Zernike length l/£ (oz) oc yjxc - X (dotted line, see Eq. (12)) 
are displayed. Within the range of values of x leading to an 
oscillatory decay, depicted by the grey regions in Fig. 5, one 
has k w l/^ oz \ Within the range of monotonic decay the 
decay rate of the leading contribution to Gij(r — > 00) is given 
by k" whereas that of the subdominant contribution is k%. 
For Xc~X < 1-3 X 1(T 3 the decay rates are fc" « l/£ (OZ) and 
k'2 ~ k, whereas for Xc — X > 4.3 x 1CP 3 the decay rates are 
fci'wKand *£'«l/£< oz >. 



a subdominant contribution oc exp(— nr-s/l — g 2 )/r (sec 
Fig. 6), i.e., as anticipated above, the leading decay at 
large distances is governed by the vicinity to the criti- 
cal point, whereas the ion-solvent coupling manifests it- 
self in the corrections to the leading behavior. Further 
away from the critical point the leading contribution de- 
cays oc exp(— k"r)/r with a subdominant contribution 
oc exp(— k'^rj/r (see Fig. 6). The inset of Fig. 6 dis- 



plays the absolute value of the real parts \k' v \ = k[ of the 
poles of Gij(k), which is identical to the wave number 
k[ (see Eq. (15)) of the oscillatory part of Gij(r) and 
which is non-zero within the grey region of Fig. 5(b). 
For the strongly antagonistic salt with / + = —20 and 
/_ = 20, in Fig. 5(b) the shortest wave length of the 
oscillations is given by (27r/fc^) m j n rj 513 (see the inset 
in Fig. 6). The corresponding value of k'{ is ps 0.0873 



so that (27r/fc^) r 



i.e. 



Gij(r) oc sm(k[r + 



const) exp(— fc"r)/r decays already within 1/45 of a pe- 
riod. In less extreme cases of solubility contrasts f± , such 
as those in Figs. 5(c) and (d), the shortest wave lengths 
are even larger. Therefore, within LDA, the oscillations 
in Gy(r), if they occur, are not expected to be experi- 
mentally detectable. In contrast, as already mentioned 
above, within BCA it is possible that the shortest wave 
lengths of the oscillations are of the order of the par- 
ticle size; such an asymptotic oscillatory decay can be 
expected to be visible in the pair distribution function. 
However, we are not aware of any experimental reports of 
Kirkwood crossover lines, which is in line with the results 
obtained within the LDA. 



IV. SEMI-INFINITE SYSTEMS BOUNDED BY 
A PLANAR WALL 

Interfacial properties such as number density profiles 
and the excess adsorption are experimentally accessible, 
albeit requiring more effort than determining bulk prop- 
erties. However, if available, they provide information 
which is significantly more detailed than the one which 
can be inferred from the bulk properties discussed in the 
previous Sec. III. In this respect the simplest and ex- 
perimentally appealing setting is that of a semi-infinite 
planar system, on which we shall focus in the follow- 
ing. Predictions within the proposed LDA (see Sec. II) 
will be compared with those obtained within the BCA as 
well as with the limited amount of corresponding exper- 
imental data which are presently available. Additional 
experimental settings are proposed in order to assess the 
predictive power of the LDA. 



A. Profiles 

According to the fluctuation-dissipation theorem the 
linear density response of a system to a weak external 
field is determined by the two-point correlation functions 
of the unperturbed system [22]. Therefore, the number 
density profiles far from a wall exhibit asymptotically the 
same type of decay towards their bulk values, i.e., either 
monotonically or damped oscillatorily, with the same de- 
cay length and periodicity as the two-point bulk corre- 
lation functions Gjj(r). (Note that this correspondence 
does no longer hold in the presence of algebraically decay- 
ing interaction potentials [39], which wc do not consider 
here.) However, from this argument one cannot draw 
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reliable conclusions concerning their structure close to 
the wall. Therefore in this latter range we determine 
the structure numerically for particular sets of parame- 
ters. Moreover, close to the wall packing effects due to 
the finite size of the fluid particles lead to layering which 
extends a few particle diameters into the system. How- 
ever, this kind of structure is not captured by the present 
square-gradient model (Eq. (1)). 

The solid lines in Fig. 7 correspond to the composition 
4>{z) [(a)], the electrostatic potential ip(z) = ip(z)kBT/e 
with e(cj))tp'(z) = -AttI b D(z) and tp(z ->• oo) -> 
[(b)], the cation number density Q+(z) [(c)], and the an- 
ion number density g~(z) [(d)] in a semi-infinite system 
bounded by a wall positioned at z = with surface charge 
density a and surface field strength h, as obtained from 
numerically minimizing the density functional in Eq. (1). 
The solvent permittivity is chosen to resemble that of 
a mixture of 3-methylpyridine (component A, ea = 10) 
and water (component B, eb — 80). The composition 
profiles 4>{z) in Fig. 7(a) turn out to be monotonic for 
weak (h = ±0.01) as well as for strong surface fields 
(h = ±1). Due to the negative surface charge density 
a, the monotonic electrostatic potential profile ip(z) in 
Fig. 7(b) is negative with surface potentials ijj(0) of some 
tens of mV, which is a common order of magnitude [40]. 
Within the range < z < 2 close to the wall the elec- 
trostatic potential ip becomes less negative upon chang- 
ing the surface field strength from h = 1 to h = — 1 
due to the increase of the permittivity as a result of the 
increase of the volume fraction 1 — of component B 
close to the wall. Similarly, due to the negative surface 
charge, close to the wall the number density g + of the 
cations in Fig. 7(c) is larger and the number density Q- 
of the anions in Fig. 7(d) is smaller than in the bulk. 
Upon changing the surface field strength from h = 1 to 
h = — 1, close to the wall the number density g + of the 
cations decreases and that of the anions, increases. 
This feature follows partly from the variation of the elec- 
trostatic potential ip. In addition, for the current choice 
of parameters the anions dissolve better in component B 
than in component A of the solvent (/_ > 0), such that 
the component B enriched near the surface (see Fig. 7(a)) 
mediates a certain preference of the anions for the wall. 

The dashed lines in Fig. 7(a) correspond to the approx- 
imate profile 



(z) = 4> b + 



C G1 



sinh((z + z )/C) ; 



GL 



sign(/i) 




(18) 



with the extrapolation length zq . Here and in the follow- 
ing we call £ = £( oz ) (see Sec. IIIB) the bulk correlation 
length. As noted in Sec. IIIB, close to the critical point 
l/£ corresponds to the asymptotic decay rate k" of the 
solvent structure factor G<p<p(r — > oo) oc exp(— r/£)/r. 
The profile 4>{z) is the analytic solution of the semi- 
infinite Ginzburg-Landau equation [41] obtained from 



minimizing f2[0, g±] after expanding Eq. (1) up to fourth 
order in — fa and neglecting the ion-solvent coupling, 
i.e., assuming f± = 0, which implies V±(0) — 0. Within 
Ginzburg-Landau theory the extrapolation length zq is 
fixed by the boundary condition 0'(O) = —3h/x- For 
\h\ — >• oo or £ — > oo this leads to \4>(0)\ — > oo. However, 
within the context of the present study 4>{z) is an ap- 
proximation of the volume fraction of solvent component 
A, which is restricted to the interval [0,1]. Hence, we 
accept the extrapolation length z as determined from 
the boundary condition 0'(O) = —3h/x only if this leads 
to (f>(0) £ [0, 1]. Otherwise the extrapolation length zq is 
inferred from the ersatz boundary condition 0(0) = 1 if 
ft > or 0(0) = if h< 0. 

The dashed line in Fig. 7(b) is the approximate elec- 
trostatic potential 



ip(z) = 4 artanh(Cps exp(— nz)), 

2ir£ B a 

£(<j) b )K 



Cpb '■= tanh ( — arsinh 



(19) 



which is the analytic solution of the semi-infinite Poisson- 
Boltzmann equation [42] for a uniform permittivity e(0fc) 
and for neglecting the ion-solvent coupling (i.e., for f± = 
0). 

Finally, the dashed lines in Figs. 7(c) and (d) are the 
approximate number density profiles of the ±-ions, 

~g ± {z) = f 6 exp(-(±^(«) + 14 (0»)) + V±{4> b )), (20) 

which correspond to the Boltzmann distributions of non- 
interacting particles in the external fields due to the ap- 
proximate electrostatic potential ip(z) and the approxi- 
mate composition <f>(z); lb = g±(z — > oo). Whereas the 
composition profile <fi(z) (Eq. (18)) and the electrostatic 
potential profile ■ip(z) (Eq. (19)) arc independent of the 
ion-solvent coupling V±, the ion number density profiles 
g±(z) (Eq. (20)) are not. 

At distances from the wall of more than a few parti- 
cle diameter (z > 2) the approximate profiles (f>(z), ^p{z), 
and ~Q±(z) differ only slightly from the ones obtained by 
the full numerical minimization. Closer to the wall the 
deviations between the numerical and the approximate 
profiles are more pronounced, but in this spatial range 
the present local model is not conclusive because it ne- 
glects the surface layering of actual fluids. A similar sit- 
uation occurs within BCA [19] (see Fig. 3(a) therein) 
shown in Fig. 8 . There, at distances z < 2A, the solvent 
composition profile 0(2) for strong ion-solvent coupling 
(/ + = 30, /_ =0, solid line) differs strongly from that in 
the absence of ion-solvent coupling (/ + = /_ =0, dashed 
line). At large distances the deviations are small. A 
closer comparison between Fig. 8 and Fig. 7(a) would re- 
quire the knowledge of the particle size, which is however 
not specified in Ref. [19]. A description of the presently 
considered semi-infinite planar system within RPA has 
been given in Ref. [9]. There the ion-solvent coupling 
has been treated perturbatively, but the full, numerically 
determined profiles 4>(z), ip(z), and g±(z) within RPA 
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FIG. 7: Profiles of the volume fraction of solvent component A [(a)], the electrostatic potential i/> (with i/>(z — > 00) — > 0) 
[(b)], the cation number density q+ [(c)], and the anion number density Q- [(d)] in a semi- infinite system bounded by a wall 
at z = z/a = with surface charge density a = — 1 ^iC/cm 2 and surface field strength /i. These results correspond to Gibbs 
free energies of transfer /+ = 0, /_ = 20, the bulk volume fraction <f>t, = 0.5 of solvent component A, and the bulk ionic 
strength / = g±b — 10 mM. The Flory-Huggins parameter x(T) is chosen to correspond to that temperature, for which the 
bulk correlation length £ is half of the Debye length 1/k, which is taken to be temperature independent (see Fig. 6). For the 
specified surface fields h the solid lines are the numeric solutions obtained from the density functional model in Eq. (1) within 
LDA. For reasons of clarity in (b) and (c) the full lines for h — ±0.01 are not designated; they can be nonetheless identified in 
an obvious way. The dashed lines correspond to the approximate profiles <p(z), tp(z), and ~Q±(z) introduced in Eqs. (18), (19), 
and (20), respectively. Note that 4>(z) and, due to the choice /+ = 0, f? + (z) are independent of the magnitude \h\; therefore 
both in (b) and (c) there is only one dashed line. For z > 2 the approximate profiles differ only slightly from the ones obtained 
by a full numerical minimization. Density oscillations close to the wall, which are expected in actual fluids, do not occur, 
because packing effects are not captured by the present square-gradient approach. 



have not been discussed. However, by comparing these 
profiles, as obtained within RPA, with those obtained 
within BCA or LDA one could assess the influence of 
non-locality on the interfacial structure in the complex 
fluids studied here. 



B. Critical adsorption 

Here we investigate critical adsorption at a wall with a 
strong surface field h and with surface charge density a. 



We consider the case that in the bulk the binary liquid 
mixture is at the critical bulk composition fa = <j) c in the 
presence of salt with bulk ionic strength g±t = I = I c 
(see Subsec. Ill A). A surface field h > (h < 0) favors 
the adsorption of A (B) particles and leads to a local seg- 
regation. In order to obtain an analytical expression for 
the excess adsorption r(£), which captures the full mean- 
field behavior to leading order close to the critical point 
(£ — s- 00), we expand the density functional in Eq. (1) in 
two steps in order to derive a Ginzburg-Landau-typc de- 
scription. In the first step the density functional ft[4>, g±] 
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FIG. 8: Composition profiles <j>(z) within BCA for strong ion- 
solvent coupling (/+ = 30, /_ = 0, solid line) and in the 
absence of ion-solvent coupling (/+ = /_ = 0, dashed line) 
taken from Fig. 3(a) in Ref. [19]. Here a solvent with ea — 80 
and eb = 20 at bulk composition cj>b = 0.09 is considered. 
The bulk ionic strength is / = 0.1 mM and the surface charge 
density is a = — 16/iC/cm 2 . At distances z < 2 A from the 
wall the two curves differ strongly from each other, whereas 
the differences are small at large distances. 



in Eq. (1) is expanded up to second order in the devi- 
ations Ag±(z) := g±(z) — I of the ion densities from 
their bulk equilibrium values I = I c . This leads to a 
density functional £li[<fi,Ag±\. Minimizing Cli[<fi,Ag±] 
with respect to Ag± renders Eulcr-Lagrangc equations 
linear in Ag±(z), the solutions of which are functionals 
Ag±(z, [<fi]) of the (up to here unknown) solvent composi- 
tion profile </>. Inserting the solutions Ag*±{z, [</>]) into the 
density functional, fii [(f), Ag±] and, as the second step, 
expanding S7i [0, Ag± [</>]] up to fourth order in the or- 
der parameter deviations ip(z) := <fi(z) — <p c leads to the 
Ginzburg-Landau-typc functional 



H[<p] 
A 



dz (a(p{z)) 2 + b( V {z)) 4 + c(<p'{z))'' 



+U{z)ip(z) -^(0) + O((A 7 )^), 



(21) 



where A is the surface area in units of a 2 . Here a = Xc~ Xi 
b = 4/3, c = x/6, and the effective "external" field is 

kctA 7 2Tre B a 2 e'((j) c ) 

U(z) = — cxp(-Kz) — exp(-2Kz). 

2 [£{<Pc)Y 

(22) 

The external field U(z) describes the influence of surface 
charges a on the order parameter tp. The first term on 
the right-hand side of Eq. (22) is due to the ion solubility 
whereas the second term is due to the dielectric proper- 
ties of the solvent. Solving perturbatively to first order 
in U the Euler-Lagrange equation, obtained from T-L /A in 
Eq. (21), leads to the equilibrium order parameter profile 

fcqiz-,0- 



From Eq. (21) one obtains the bulk correlation length 
£ = A + 0((A 7 ) 2 ) = X^~7 + G((A 7 ) 2 ) = e oz) + 



0((A 7 ) 2 ) with A = 



X 



(see Eqs. (9), 



6(Xc - X) V « 
(10), and (12) with Xs{4>cJc) = Xc) for X < Xc- In 
the following we shall neglect the corrections 0((A 7 ) 2 ), 
which are expected to be small within LDA (see Sec. III). 



Using the empirical form x(T) = xs 
one obtains £(T -> T c ) '4° $\t\~ l 



XH 

T 



(see Sec. II) 



XcTc 



with the non- 
the critical 



universal critical amplitude £j~ = , 

(mean-field) exponent v = 1/2, and t = (T - T c )/T c 
Moreover, from Eq. (21) one obtains the bulk order pa 



rameter tp e q{z 



»;0 = \l~2b = 



3(x~Xc) t-^o 



8 



to 



t\P with the critical amplitude 



to 



isIxhI 



and 



)j 8T C 

the critical (mean- field) exponent (3 = 1/2. For later 
purposes (see the text below Eqs. (23) and (26)) here 
we note that \fcjb = \J2xd^ = V^mo^o =1/2 for 
X = Xc = 2 (see Sec. Ill A). 

At the critical point (£ = oo) and far away from 
the substrate the equilibrium order parameter profile 
<^ oq (z;£ = oo) decays as 



oo; £ = oo) : 
sign(/i) 



sign(fe) 

2z 
9ctA 7 
' 10k 3 " 



(23) 



9n£ B a 2 e'(<f) c )\ 1 



4O K 4 (e(0 c )) 2 ) z 



+0(z~ 3 ). 

The leading contribution sign(/i) — 

the scaling form sign(h)moc+(z/^Q ) ~&' v with the univer- 
sal amplitude c + = y/2 [43] , where the critical exponents 



can be written in 



take their mean-field values f3 



\ [33]. Accordingly, 



the leading term in Eq. (23) is not affected by the sur- 
face charge, the presence of ions, or the dielectric prop- 
erties of the solvent. However, these materials properties 
do modify the amplitude of the subleading contribution 
(oc 1/z 2 ). 

Close to the critical point (£ — > oo) the excess ad- 

oo 

sorption r(£) = J dz (p cq (z;£) with the perturbatively 
o 

obtained profile ip C q(z',£) (see above) is given by 



r(0=r„(0 + ri(0 + 0(1/0 



with 



and 



r (0 



sign(/i) 



ln(20 



ri(0 



3crA 7 3ire B a 2 e'(<l) c ) 



8k 2 



(24) 



(25) 



(26) 
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FIG. 9: Comparison of the numerically calculated excess ad- 
sorption r(£) obtained within the full model in Eq. (1) (•) 
with the predictions of Eq. (24) for the parameters used in 
Fig. 7 with h = 1. The Debye length k _1 (marked by an 
arrow) corresponds to a bulk ionic strength I = 10 mM. The 
term To(£) (dotted line, see Eq. (25)), which contains the 
leading contribution to the excess adsorption and which corre- 
sponds to a vanishing surface charge density (a — 0), exhibits 
visible deviations from the numerical results (•); nonetheless 
Fi(£)/ro(£) — > for £ — > oo. Taking into account in addition 
the term Ti(^) (see Eq. (26)), which exhibits a dependence on 
the surface charge a, quantitative agreement is found between 
To(5) + Ii(£) (solid line) and the numerical results (•) in the 
limit £ — > oo. This finding also implies that those terms of 
Eq. (1), which have been left out upon deriving Eq. (21), do 
not contribute detectably. 



The leading contribution sign(/i)W - ln(£) can be written 

in the scaling form ~ sign(/i)mo£o <7+(— ln(|t|)) with the 
universal amplitude g + = \[2v = l/\/2 within mean- 
field theory [43]. To(^) diverges for £ — > oo whereas 
Ti(^) remains finite and thus represents the first sub- 
dominant correction. Figure 9 compares the predictions 
of Eq. (24) with the results obtained by numerically cal- 
culating the excess adsorption within the full model as 
given by Eq. (1) (•) for the parameters used in Fig. 7 
with h = 1, in particular for large £. Whereas the leading 
contribution ro(£) in Eq. (25) itself (dotted line) deviates 
visibly from the full numerical results (•), there is quan- 
titative agreement between the latter and Tq(£) + Ti(£) 
(solid line) in the limit £ — > oo. Since To(^) corresponds 
to a vanishing surface charge (er = 0), the difference be- 
tween the dotted and the full line in Fig. 9 demonstrates 
the influence of electrostatic interactions on the excess 
adsorption. The quantitative agreement of Tq(£) + Ti(^) 
(solid line) with the numerical results (•) indicates that 
the terms of Eq. (1) neglected upon deriving Eq. (21) 
do not contribute detectably to the leading and the first 
sublcading behavior of r(£ — s- oo). Moreover, for the 
given choice of parameters the magnitude of the correc- 



tion 1 0( 1/£)| in Eq. (24) turns out to be smaller than 
|ri(£)|, which in turn vanishes relative to T(£ — > oo). 
Note that within BCA, for /+ = 0, /_ = 20 (as used 
in Fig. 9) the uniform bulk state is unstable in a cer- 
tain vicinity of the critical point, because \g\ = 1.39 > 1 
(see Sec. IIIB), which precludes calculating the excess 
adsorption T(£). 

Critical adsorption occurs upon approaching the criti- 
cal point (£ — > oo), where the excess adsorption diverges 
as r oc ln£, which is in agreement with the expected 
universal scaling behavior V oc £ 1_ < 3 / l/ [43, 44] for the 
classical exponents f3 = v = 1/2 corresponding to the 
present mean-field theory. It is apparent from Eq. (24) 
that the leading contribution r is not altered by the sur- 
face charge, the presence of ions, or the dielectric prop- 
erties of the solvent. However, these non-universal prop- 
erties do influence the subleading contribution Y\. 

Recently the adsorption of critical water+2,6- 
dimethylpyridine mixtures with KBr of various ionic 
strengths / has been investigated by means of surface 
plasmon resonance [34]. For the case of a hydrophobic 
wall the excess adsorption turned out to be practically 
independent of the ionic strength. This is in agreement 
with Eq. (24) because a hydrophobic wall is only weakly 
charged [45] such that the second and third terms on the 
right-hand side of Eq. (24) are negligibly small. 

For the case of a hydrophilic, negatively charged (a < 
0) wall a decrease of the adsorption of water has been 
measured upon adding salt [34]. Hydrophilic walls can 
be expected to be strongly charged [46] such that a = 
sign(cr)f7 sat with the saturation surface charge density 
cr sat = K£(0 c )/(7r£ B ) [47]. In this case from Eq. (24) 
one obtains 

9T 3 / . . , A s'(4> c )\ 
dI=-^H^ + l^)- (27) 

Equation (27) assumes only a weak dependence of cf> c 
(and thus of A7 and of e(4> c )) on the ionic strength I (see 

Sec. Ill A) so that the derivative does not appear. 

al 

If 2,6-dimcthylpyridine is denoted as the A component 
and water as the B component of the binary liquid mix- 
ture (i.e., r measures the excess of 2,6-dimcthylpyridinc), 
at the lower critical demixing point an experimental value 
of e'/e ps —1.2 is found [48, 49]. For this mixture the sol- 
ubility contrasts for KBr are /+ ps 2.5 and /_ ps 8.4 [24] 
which leads to /+ - /_ ps -5.9 and A7 ps -0.30. Within 
LDA from these numbers one finds dT/dl > (i.e., de- 
creasing water adsorption upon adding salt), and the sec- 
ond (dielectric) contribution on the right-hand side of 
Eq. (27) dominates. In contrast, within BCA one has 
<9p(BCA)^j < (i.e., increasing water adsorption upon 
adding salt), because A7 ps / + — /_ leads to a dom- 
inance of the first (ion solubility) contribution on the 
right-hand side of Eq. (27). Hence the overestimation of 
the ion-solvent coupling within BCA leads to a sign of 
c?p(BCA) jdl which is not compatible with the aforemen- 
tioned experimental findings in Rcf. [34] , whereas the sign 
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of dT/dl within the present LDA is in agreement with 
these findings. Since the dielectric properties are exper- 
imentally accessible one could use Eq. (27) to determine 
A7 from measurements of the excess adsorption F as a 
function of the ionic strength /. A comparison of this 
resulting value for A7 with the difference /+ — /_ of the 
Gibbs free energies of transfer (inferred, e.g., from elec- 
trochemical methods) would be a direct way to probe 
quantitatively the difference between the LDA and the 
BCA. 

In order to further test the different predictions fol- 
lowing from BCA and LDA we suggest additional ad- 
sorption measurements for hydrophilic walls. For KBr 
as salt the arguments above lead to the assertions that, 
within LDA, one has dY/dl > (i.e., decreasing water 
adsorption upon adding salt) independently of the sign of 
the surface charge a (because the last term on the right- 
hand side of Eq. (27) dominates), whereas within BCA 
(9r( BCA Vd/ is expected to change sign upon changing 
the sign of a (because the first term on the right-hand 
side of Eq. (27) dominates). More interestingly, using 
an antagonistic electrolyte (i.e., with /+ and /_ having 
opposite signs) such as HBr (/+ m -11.2, /_ m 8.4 [24], 
i.e., A7 rj —4) the first (ion solubility) contribution on 
the right-hand side of Eq. (27) is dominating such that 
dT/dl and a are expected to have the same sign. In 
this case, upon adding salt, the amount of adsorbed wa- 
ter either decreases or increases depending on the sign of 
the surface charge. Fhis is in contrast to electrowetting 
where the water adsorption increases with the magnitude 
but independent of the sign of the surface charge [50]. 
Within the BCA approach of Ref. [21] the difference for 
cations and anions with respect to their solubility con- 
trasts in the two pure solvent components is neglected 
(i.e., /+ = /_ >• 1 so that A7 = 0) to the effect that the 
reported capillary condensation-like adsorption of water 
between two equally charged walls at variable distance 
should be independent of the sign of the surface charge. 
The analysis above implies that the same property is ex- 
pected to occur for sufficiently small values of IA7I, but 
the adsorption may depend on the sign of the surface 
charge if | A7I becomes of the order unity. 



V. COLLOIDAL INTERACTIONS IN 
NEAR-CRITICAL ELECTROLYTE SOLUTIONS 

The interactions between colloidal particles in a fluid 
medium comprise dispersion forces, direct screened 
Coulomb forces, steric forces, as well as solvent-mediated 
interactions, which are commonly referred to as solvation 
forces. If the thermodynamic state of the fluid medium is 
moved towards a critical point, e.g., the critical demixing 
point of a binary liquid mixture, the fluctuation induced 
long-ranged, and universal critical Casimir force emerges; 
this singular contribution to the solvation force domi- 
nates [51-53]. In the presence of a sufficiently strong ad- 
sorption preference of colloids or of confining walls for one 



of the components of the solvent the critical Casimir force 
depends only on the relative signs of the surface fields 
h acting on the order parameter at these surfaces (and 
on the geometry of the latter) but not on non-universal 
material parameters. The critical Casimir force is attrac- 
tive (repulsive) if the surface fields have equal (opposite) 
signs. 

While being investigated theoretically for quite some 
time [51], the direct experimental verification of the crit- 
ical Casimir effect has been achieved only recently for a 
single colloidal particle dispersed in a mixture of water 
and 2,6-dimcthylpyridine close to a planar wall [53-55]. 
In that study the measured effective colloid- wall interac- 
tion potential was interpreted in terms of a superposi- 
tion of a repulsive screened Coulomb force and the criti- 
cal Casimir force. Within this picture the direct electro- 
static repulsion dominates for temperatures T which de- 
viate from T c more than a few tenth of a Kelvin, whereas 
an increasingly strong Casimir attraction (repulsion) oc- 
curs for symmetric (antisymmetric) boundary conditions 
upon approaching the critical point (T — > T c ). Recently 
these experiments have been modeled within RPA, which 
provides a satisfactory fitting of the experimental curves 
in Refs. [54, 55]. However, this approach involves a large 
number of model parameters [56], which limits the con- 
clusiveness of these fits. 

The influence of adding salt onto the effective colloid- 
wall interaction has been studied recently [34] using the 
same experimental setup as in Refs. [54, 55]. It turns 
out that with 10 mM of KBr the Casimir attraction for 
symmetric boundary conditions starts to dominate the 
direct electrostatic repulsion already several Kelvin away 
from the critical point (instead of tenths of a Kelvin as 
for the salt-free solvent). Moreover, for antisymmetric 
boundary conditions, for which both the direct electro- 
static and the critical Casimir forces are expected to be 
repulsive, in the presence of salt an attraction has been 
detected within an intermediate temperature range. The 
latter observation demonstrates that, under certain con- 
ditions, assuming the simple superposition of the direct 
electrostatic and the critical Casimir forces is insufficient 
to understand the actual effective interaction. 

Based on the Ginzburg-Landau-likc description in 
Eq. (21), which follows from the full model given in 
Eq. (1), we have identified a mechanism giving rise to the 
aforementioned unexpected attraction for antisymmetric 
boundary conditions in the presence of salt [26]. For 
A7 7^ the cations and anions are separated close to the 
surfaces, where the order parameter tp eq is non- uniform. 
This gives rise to dipolar layers, which can interact with 
the surface charges at the distant surface. These dipolc 
layers are expected to contribute significantly to the ef- 
fective interaction if the direct electrostatic interaction is 
weak; this is the case for antisymmetric boundary con- 
ditions, for which the hydrophobic wall is expected to 
be weakly charged. If the direct electrostatic interac- 
tion is strong, which is expected to occur for symmetric 
boundary conditions with hydrophilic surfaces, the dipo- 
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lar layers do not significantly contribute to the effective 
interaction. In that case the sole effect of adding salt is 
to reduce the Debye length and thereby to weaken the di- 
rect electrostatic repulsion relative to the critical Casimir 
force. Therefore, with salt the onset of effective attrac- 
tion occurs at temperatures T further away from T c than 
in the salt- free case [26] . 

A conceivable alternative mechanism for the emergence 
of an effective attraction in the case of antisymmetric 
boundary conditions has been proposed which is inde- 
pendent of differences in the solubility of cations and an- 
ions [57]. It has been argued within RPA that a charged 
wall (of either polarity) accumulates an increased num- 
ber of ions compared to an uncharged wall. Due to this 
enhanced total density of ions (which are hydrophilic in- 
dependent of their sign) a charged wall should, from a 
distance, appear increasingly hydrophilic upon adding 
salt such that for certain system parameters an under- 
lying actually hydrophobic character of a wall might be 
overcompensated and turn into an effectively hydrophilic 
wall; actually hydrophilic walls remain so upon adding 
salt [57]. Such a salt-induced apparent hydrophilicity, 
which would occur on the surfaces of all dissolved col- 
loids, would in turn lead to effectively symmetric bound- 
ary conditions and thus to attractive solvation forces. 
However, according to Fig. 9, even in the presence of salt 
the excess adsorption follows the actual preference of the 
surface field, also upon approaching T c . As mentioned 
in Sec. IV B, it is not possible to calculate the excess 
adsorption within BCA (and RPA) for the parameters 
used in Fig. 9, because for them within this approach 
the uniform bulk state is thermodynamically unstable. 
Thus Fig. 9 demonstrates that within LDA salt-induced 
apparent hydrophilicity does not occur (i.e., T does not 
become negative). Therefore there is reason to expect 
that salt-induced apparent hydrophilicity is an artifact 
of the BCA and the RPA. In addition, by means of sur- 
face plasmon resonance it has been checked experimen- 
tally that the adsorption preference, in particularly of hy- 
drophobic substrates, is not altered by adding salt [34]. 
Therefore, for antisymmetric boundary conditions of the 
order parameter, there are doubts that salt-induced ap- 
parent hydrophilicity can serve as an explanation for the 
experimentally observed effective attraction within an in- 
termediate temperature range. 

Recently, additional numerical studies within BCA 
have been performed suggesting that ion- induced "pre- 
cipitation" [8] or non-linearities [58] influence the ef- 
fective colloid-colloid interaction. However, our results 
in Sees. Ill and IV concerning the reliability of BCA 
point towards the possibility that those proposed effects 
are artifacts of the BCA due to an overestimation of 
the ion-solvent coupling. Therefore it would be worth- 
while to reconsider the aforementioned proposed effects 
within LDA, which has been shown to be consistent with 
presently available experimental evidence. 



VI. CONCLUSIONS AND SUMMARY 

We have derived a local density approximation (LDA) 
for the density functional of point-like ions interacting 
locally with a binary liquid mixture acting as a solvent 
(Fig. 2), which for large free energies of transfer of the 
ions improves the frequently used bilinear coupling ap- 
proximation (BCA). It turns out that within the pro- 
posed LDA, and in contrast to the BCA, the influence of 
ions on the bulk phase diagram (Fig. 1), the critical point 
(Fig. 3), and the bulk structure (Figs. 4-6) is predicted 
to be weak. This is in agreement with the presently avail- 
able experimental data. The interfacial structure at dis- 
tances from a charged wall less than a few particle diam- 
eters can be predicted by neither BCA nor LDA, because 
none of these two local models accounts for the layering 
due to packing effects which dominate in actual fluids at 
such short distances. But, within both BCA and LDA, 
further away from the wall the system can be described 
reliably in simple terms of a uniform permittivity and 
an effective surface field. This description nonetheless 
captures the actual structure as obtained from the full 
numerical minimization (Fig. 7). Upon approaching the 
critical point the subleading (but not the leading) con- 
tribution to critical adsorption is found to be sensitive to 
system and materials parameters such as the bulk ionic 
strength, solubility properties, surface charges, and the 
permittivities of the solvent components (Fig. 9). 

If a salt disturbs the molecular arrangement of the sol- 
vent molecules only at small distances, it is modeled here 
by point-like ions which interact locally with the solvent. 
Such a salt does neither significantly alter the bulk phase 
behavior nor the bulk structure or the asymptotic decay 
of density profiles at walls. But it contributes to the in- 
terfacial structures up to distances of the order of the 
Debye length as well as to critical adsorption in sub- 
leading order. If the solvent is moved thermodynamically 
towards a critical point, where the bulk correlation length 
£ diverges, the ratio l/(re£) of the range 1/re of the ion- 
related surface structure and £ becomes small, such that 
the leading, universal critical behavior of the solvent is 
not altered by adding salt [26]. 

Whether the description of a given electrolyte solution 
within a local model is justified or not does not depend 
on the salt alone but on the combination of salt and 
solvent. Experimentally observed effects in binary liq- 
uid mixtures due to adding salt, such as the shift of the 
critical point, depend sensitively on the type of mixture 
(compare Rcf. [31] for water+2,6-dimethylpyridine and 
Ref. [14] for heavy water+3-methylpyridine). Moreover, 
the measured critical point shifts exhibit a strong depen- 
dence on the size of the ions (compare Ref. [14] for al- 
kali halides and Ref. [15] for sodium tctraphcnylborate). 
These evidences in combination with the present analysis, 
which implies a weak influence of the ionic charge, lead to 
the conclusion that stcric effects might play an important 
role for the ion-solvent interaction. This interpretation 
is supported by reports of critical point shifts of similar 
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magnitude in binary liquid mixtures due to adding non- 
ionic impurities [59]. Consequently it appears as if it is 
mostly the property of an ion to be a structure maker or 
a structure breaker and only to a lesser extent its elec- 
tric charge which determines the influence of a salt onto 
the properties of the solvent. In order to obtain quan- 
titatively reliable predictions for electrolyte solutions, a 
larger effort has to be devoted to study the steric and 
chemical influence of ions on the solvent. 

In summary, for ions dissolved in a binary solvent 
the bulk phase diagrams (Fig. 1), the critical point 
shifts (Fig. 3), the bulk two-point correlation functions 
(Figs. 4-6), the surface structures (Fig. 7), and critical 
adsorption (Fig. 9) have been investigated within a novel 
local density approximation (LDA) for the ion-solvent 
interaction (Fig. 2). The commonly used bilinear cou- 
pling approximation (BCA) turns out to strongly over- 
estimate the influence of ions on the solvent properties 
in cases of realistic solvation free energies, whereas the 
LDA introduced here, in agreement with various avail- 
able experimental data, predicts small effects. Although 
the presented LDA is expected to be more accurate than 
the BCA, the former requires the same and not more pa- 
rameters than the latter. According to its derivation, the 
LDA is expected to be reliable for small ionic strengths 
and on length scales larger than the particle size. Both 
available and possible future experiments have been dis- 
cussed to probe and to explore the reliability and the 
range of validity of the present theoretical approach. 
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Appendix A: Derivation of the solvent-induced ion 
potential V±{cj>) 

Within the present model the solvation of ions is dom- 
inated by short-ranged interactions with solvent parti- 
cles. Accordingly, the ion-solvent interaction is mod- 
eled locally in terms of a free energy density QiVi(4>) 

i=± 

(see Eq. (1)). In order to derive an expression for the 
solvent-induced ion potential V± (</>), resorting to a lat- 
tice gas model and in line with the short range of the 
ion-solvent interaction, a single site is considered which 
is occupied by one solvent particle (either of type A or 
of type B, i.e., corresponding to occupation numbers 
Na = 1 — N B € {0, 1}) and an arbitrary number of pos- 
itive and negative ions (N± S No). The surrounding of 
this considered site acts as a particle reservoir which is de- 
scribed by chemical potentials Xik B T,i € {A, £?,+,— }, 
and the interaction energy per k B T of two particles of 



species i,j € {A,B,+,— } on that site is given by kij 
with = for i,j 6 {A, B}, because a particle does 
not interact with itself and the site cannot be occu- 
pied by more than one particle of type A or B. The 
chemical potentials A^ are related to the chemical po- 
tentials /i0 and fi± introduced in Sec. II (see below). 
Accordingly, the Hamiltonian Hk B T of a configuration 
(N A ,N B ,N+,N-) G {0,1} x {0, 1} xN xN on one site 
with Na + Nb = 1 is given by 

H(N A ,N B ,N+,N-) = k A+ NAN + + kA-N A N_ + 

k B+ N B N + + k B -N B N- + 
H ion (N+,N_) (Al) 

with the ionic part 

H ion (N + ,N_) = (A2) 

, 2V + (iV + -l) , Ar Ar , N-(N--1) 
k++^^ '- + k + _N + N_ + k__ K — '-. 

The corresponding grand partition function is 

C(a a ,a b ,a + ,a_)= J2 ( A3 ) 

(N A ,N B ,N + ,N_) 

exp( J2 \iNi-H(N A ,N B ,N + ,N-j), 

ie{A,B,+,-} 

where the outermost summation is over all configurations 
(N A ,N B ,N+,N_) G {0,l}x{0,l}xN o xN o with N A + 
N B = 1. With N A =Nt,N B = l- N^ for N^ G {0, 1} 
the Hamiltonian in Eq. (Al) can be rewritten as 

H{N4,,1-N^N + ,N-) = f + N<t,N + + f-N^N- + 

k B+ N+ + k B _N- + 
H ion (N +: N_) (A4) 

with the solvation energy difference f± = kA± — k B ± 
of a ±-ion. The grand partition function in Eq. (A3) is 
C(Aa,A b ,A + , A_) = exp(A B )Z(/i0,^ + ,^_) with 

^W)M+>M-)= X! exp ( X] 

- f+N^N + - f-N^N_ - H ion (N + , N-)J , (A5) 

where ^ = A^ — X B and fi± = \± — k B ± (see Sec. II) and 
where the outermost summation is over all configurations 
(JV ,JV + ,AT_)g{O,1}xNoxNo. 

In the limit of low ionic strength (/z± — > — oo; this 
implies A± — > — oo which in turn means that the averages 
(N±) = g± are small) one obtains 

~ 1 + exp(^) + (A6) 
exp(/x+)(l + exp(^ - /+)) + 
cxp(^_)(l + exp(^ - /_)). 

d In Z 

Using <t> = (N4,) = and q± = (N±) = 

^- one finds from the grand canonical potential 
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— fcs^ln Z(fx^, fi-) the Helmholtz free energy per with M± = (1 — 4>(1 — cxp(— /±))) 1 . From Eq. (A7) 
kgT one infers the effective ion-solvent interaction V±(<f>) := 



-hlZ(lJ, ( f„fX + ,fJ > -) + + Q+l-i+ + Q-fl- 

01n^+(l-0)ln(l-0)+ (A7) 
e+(ln £>+-! + In M+) + £>_ (In q_ - 1 + In M+) 
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